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ABSTRACT 

HD 15137 is an intriguing runaway O-type binary system that offers a rare opportunity to 
explore the mechanism by which it was ejected from the open cluster of its birth. Here we present 
recent blue optical spectra of HD 15137 and derive a new orbital solution for the spectroscopic 
binary and physical parameters of the O star primary. We also present the first XMM-Newton 
observations of the system. Fits of the EPIC spectra indicate soft, thermal X-ray emission 
consistent with an isolated O star. Upper limits on the undetected hard X-ray emission place 
limits on the emission from a proposed compact companion in the system, and we rule out a 
quiescent neutron star in the propellor regime or a weakly accreting neutron star. An unevolved 
secondary companion is also not detected in our optical spectra of the binary, and it is difficult to 
conclude that a gravitational interaction could have ejected this runaway binary with a low mass 
optical star. HD 15137 may contain an elusive neutron star in the ejector regime or a quiescent 
black hole with conditions unfavorable for accretion at the time of our observations. 

Subject headings: stars: individual(HD 15137), binaries: spectroscopic, X-rays: binaries 



1. Introduction 

While most 0- and B-type stars are believed 
to form in open clusters and stellar associations, 
some are observed at high galactic latitudes and 



none have been found (Philp ct al. 111996 Sayer et 



with large peculiar space velocities (Gies & Bolton 



1986). These stars were likely ejected from the 
clusters of their birth, and there are two accepted 
mechanisms to explain the origin of their run- 
away velocities. Close multi-body interactions 
in a dense cluster environment may cause one 
or more stars to be scattered out of the region 
( Poveda et al.|1967 l, and about 10% are expected 
to be ejected as binary pairs ( |Leonard fc Duncan 
1990). An alternative mechanism involves a su- 



pernova explosion within a close binary, ejecting 
the secondary due to the conservation of momen- 
tum (Zwicky 19571. The resulting neutron star 



(NS) or black hole (BH) may remain bound to the 
secondary if not enough mass is lost during the 
explosion or if a retrograde kick velocity occurs, 
and a binary fraction of 20 — 40% is predicted for 



these runaways (Portegies Zwart 2000|. The ob- 
served fraction of binaries among runaways seems 
consistent with either scenario (5 — 26% 



Mason 



et al. 19981. While both processes do occur in 
nature, it is not clear which process is dominant. 

Distinguishing between the dynamical ejection 
and binary supernova scenarios for an isolated run- 
away star can be nearly impossible. However, bi- 
nary runaway systems offer a significant advantage 
to cluster ejection studies since multiwavclcngth 
observations can pinpoint NS companions pro- 
duced in supernovae. Optical spectroscopic stud- 
ies alone offer circumstantial evidence for com- 
pact companions, primarily by high eccentricity 
and low mass function determined from the orbit, 
but also due to nitrogen and/or helium enrichment 
and fast rotational velocity that suggest prior mass 
transfer ( Hoogerwerf et al7||2001 1. A few runaway 
massive stars are predicted to be associated with 



radio pulsars (Portegies Zwart 20001, but so far 
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al. 19961. X-ray emission is a better diagnostic 
since stellar wind interactions with a NS produce 
a measurable X-ray excess above the normal, in- 
trinsic emission from shocks in the stellar winds 



(Sana et al.|2006 


Popov et al.|2000 


Lamers et al. 


1976 


)• 



HD 15137 is a runaway binary that offers the 
rare opportunity to diagnose its cluster ejection 



mechanism. Its spectral type is 09.5 III(n). Boy 



ajian et al. ( 2005 1 found that this runaway single- 
line spectroscopic binary (SB1) system was likely 
ejected from the open cluster NGC 654, but the 
travel time since its ejection, 10 Myr, presents an 
unusual paradox since it is longer than the ex- 
pected lifetime of an O-type star. |Boyajian et al 



(2005} and |McSwain et al.| fl2007a| ) presented pre 
liminary orbital solutions for the SB1, finding a 
highly eccentric orbit with a period of about 30 
d. The O-type star has an effective temperature 
T e ff = 29700 K and surface gravity logg = 3.50, is 
located at a distance of 2.420 kpc, and has a fast 
peculiar space velocity V^ cc = 63 km s _1 ( |Mc-| 



pwain et al. 2007a I. Such a value of V poc places 



this SB1 system among the class of runaway O 
stars. 

The very low mass function suggests a low mass 
companion in HD 15137, yet so far Boyajian ct al. 
( |2005[ ) andjMcSwain et al.| (|2007b|) have tried un- 



successfully to identify the nature of the compan- 
ion star. Using doppler tomographic separation 
with their preliminary orbital solution, Boyajian 



et al. ( 2005 ) were unable to identify an uncvolved 
secondary companion in their red optical spectra. 
Furthermore, we simulated 4-body interactions of 
HD 15137 with the STARLAB package and found 
a very low probability that dynamical interactions 
could have ejected a low mass, unevolved star with 
the O-type primary (McSwain et al. 2007b[). Its 
relatively short period, eccentric orbit led Boya-| 



jian et al. ( 2005 1 to suggest that the orbit was 
widened during a supernova event and that tidal 
circularization has not yet occurred. The system 
is not a known X-ray source, so it has been pro- 
posed to be a high mass X-ray binary (HMXB) 
candidate. Current upper limits of the X-ray flux 
rule out a stellar wind-accreting NS in the system, 
and we performed a radio pulsar search with null 
results ( |McSwain et~aL||2007b I . 

In this work, we present new blue optical spec- 
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tra of HD 15137 in §2. Based on new radial veloc- 
ities of the system, we argue that the orbital pe- 
riod is somewhat longer than we have previously 
determined, and we present an improved orbital 
solution for this SB1. In §3, we present the stellar 
physical parameters determined by comparing our 
optical spectra to model spectra. We describe the 
first X-ray detection of HD 15137 with the XMM- 
Newton observatory in §4. From fits of the soft 
X-ray spectra, we place upper limits on the hard 
X-ray flux from a potential NS companion. We 
discuss these results further in §5, and summarize 
the key results in §6. 

2. Optical Observations and Radial Veloc- 
ity Measurements 

We obtained 44 spectra of HD 15137 at the Ob- 
servatoire de Haute-Provence (OHP) during sev- 
eral observing runs from 2005 October to 2007 
November, including 23 spectra obtained over 30 
consecutive nights in fall 2007. We used the 1.52m 
telescope with the Aurelie spectrograph with grat- 
ing #3 and the 2048 x 1024 CCD EEV 42-20#3 
detector with a pixel size of 13.5 /im. These spec- 
tra cover a wavelength range between 4460-4890 
A with a resolving power of R — A/AA « 8000. 
The exposure times were typically 30 minutes but 
sometimes longer depending on weather condi- 
tions. On 2006 October 27, we performed an in- 
tense stare at HD 15137 with 12 consecutive expo- 
sures, although on most other nights only a single 
spectrum was obtained. The Aurelie CCD data 
consist of five images of the object's spectrum ob- 
tained using a Bowcn-Walvarcn image sheer (de- 
scribed in |Gillet et al.||1994| ). The OHP data 
were bias corrected, flat fielded, cleaned for cos- 
mic rays, and extracted using the MIDAS software 
developed at the European Southern Observatory. 
We performed wavelength calibration of each spec- 
trum using a ThAr comparison spectrum obtained 
shortly before or after each observation. The spec- 
tra were rectified using line-free continuum regions 
and interpolated onto a log wavelength scale using 
a common heliocentric wavelength grid. 

We also obtained 47 spectra of HD 15137 at 
the KPNO coude feed (CF) telescope over 35 con- 
secutive nights during 2008 October and Novem- 
ber. We used grating B in third order with 
the 4-96 order-sorting filter and the F3KB detec- 



tor. This instrumental configuration resulted in a 
wavelength range of 4130-4570 A with R w 9000 
across the chip. We generally obtained one spec- 
trum of 15 minutes duration each night, except 
for one night (2008 November 1) when we per- 
formed another long stare at HD 15137, collect- 
ing 16 spectra over four hours. We also obtained 
ThAr comparison spectra about every 2-3 hours 
throughout each night. The spectra were bias cor- 
rected, flat fielded, and wavelength calibrated us- 
ing standard procedures in IRAB^] Finally, the 
coude feed spectra were rectified to a unit contin- 
uum using line-free regions and interpolated onto 
a log wavelength scale using a common heliocen- 
tric wavelength grid. 

The recent optical spectra of HD 15137 show 
evidence of low amplitude, high frequency line pro- 
file variations, especially during the long stares 
that were performed in each of the runs described 
above. The variations in the He I A4471 line are 
shown in Figure [l] While the bumps in the line 
profile are consistent with the noise level of the 
continuum, the line profile is also non-symmetric 
and shifting slightly over the long stare. These 
short-term variations, likely due to nonradial pul- 
sations, were suggested originally by |Boyajian et 
al. (20051 and have been described in more de- 
tail by De Becker et al. (20081. Superimposed on 



these low amplitude line profile variations, we ob- 
served larger radial velocity, V r , variations of each 
absorption line. Thus HD 15137 appears to be 
a single-line spectroscopic binary (SB1) that also 
exhibits line profile variations. 

To measure V r of HD 15137, we used a cross- 
correlation procedure similar to that described in 



McSwain et al. (2007a|. However, we used a mean 



spectrum created from each data set, rather than 
any single observation, as a reference spectrum for 
the cross correlation process. This choice allowed 
us to minimize the effects of small-scale line pro- 
file variations and identify the bulk V r shifts in 
the data instead. Using this mean spectrum, we 
fit the core of each absorption line with a parabola 
to determine its absolute radial velocity. The re- 
maining spectra were then cross correlated with 
this template to determine V r for each line. We 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Sciences Foundation 
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used all available strong, unblended absorption 
lines to measure V r , and rest wavelengths for each 
line were taken from the NIST Atomic Spectra 
Databas^] In the CF data set, the lines used to 
measure V r were H7, He I AA4144, 4388, 4471, and 
He II A4200. In the OHP data set, we used the 
H/3, He I AA4471,4713, and He II AA4542,4686 
absorption lines. 

The resulting mean V r and standard deviation, 
a, are listed in Table [l] In many cases, the mea- 
sured a probably underestimate the true error 
in V r due to the line profile variations. During 
both of our long stares at HD 15137, we observed 
changes in V r of about 10 km s _1 over only a few 
hours. We recommend adding this error of ±5 
km s _1 in quadrature with a for a more represen- 
tative error in V r . 

In our previous spectroscopic studies of HD 
15137 (Boyajian et al. 2005 McSwain et al 



2007a I, we proposed an orbital period P 
However 



30 



d for the system. However, our new V r measure- 
ments exclude such a value of P, especially using 
our data from two observing runs that each took 
place over ~ 30 consecutive nights. We performed 
a new period search on all available V r measure- 
ments of HD 15137 using a version of the discrete 
Fourier transform and CLEAN deconvolution al- 



of R 



gorithm of Roberts, Lehar, & Dreher (1987) (writ 



ten in ID IP] by A. W. Fullerton). In addition, we 
used the phase dispersion minimization (PDM) 
algorithm (Stellingwer f|1978| ) to compare with the 
CLEANed power spectrum. The PDM method is 
better suited for eccentric orbits that are strongly 
non-sinusoidal. 

There was no one clear signal that stands out 
from any of the resulting periodograms, so we in- 
spected each candidate frequency carefully. We 
used each proposed period as input into the non- 
linear, least-squares fitting program of |Morbey fc| 



Brosterhus ( 1974 1 to solve for the resulting orbital 



elements. After ruling out all resulting V r curves 
with poor fits and extremely large scatter (and 
also any P < 35 d), we settled upon two accept- 
able solutions with nearly equally good fits and 
frequencies 0.018 d _1 and 0.015 d" 1 , correspond- 
ing to P = 55.4 d and P = 65.2 d, respectively. 



3 The NIST Atomic Spectra Database is available online at 



Since our data are often spaced at intervals on 
the order of 30 d, this produces a sampling fre- 
quency f s — 0.033 d _1 , which is also a signal 
found by CLEAN (and also corresponds to our ear- 
lier proposed period). Our sampling rate induces 
an alias frequency f a = f s — /, where / is the 
true signal frequency. Since no significant alias- 
ing would occur if the Nyquist condition is satis- 
fied, / < f s /2, we attribute the lower frequency 
as the alias and adopt the true orbital frequency 
/ = 0.018 d _1 . Thus we adopt an orbital period 
P = 55.4 d for HD 15137. We note that a long- 
period aliasing effect also appears prominently in 
our results; the strong peaks near 0.0015 d^ 1 and 
0.007 d" 1 found by CLEAN are the result. 

As further verification of our adopted orbital 
period, we tried manipulating our collection of V r 
data in several ways to test whether the period 
search results remained consistent. We performed 
a period search with the CF dataset omitted, find- 
ing two potential periods of 57 d and 67 d. When 
we instead omitted the OHP data from the col- 
lection, a period of 52 d stands alone as the most 
prominent signal. Finally, we repeated the period 
search an additional ten times, randomly removing 
20% of all of the V r points in each trial. Promi- 
nent signals remain consistently at 55 d and 65 
d. These results provide further support that the 
timescales of 55 d and 65 d are good candidates 
for the orbital period, with the former preferred. 
However, we recommend that both values be con- 
sidered with caution since low frequency trends 
also appear to contribute to the temporal behavior 
of the radial velocities. The corresponding orbital 
parameters should be treated with caution as well. 

Finally, we inspected the complete V r dataset 
for systematic velocity differences by repeating the 
orbital fit for each set separately, allowing only Vb 
to vary. The V r offset was determined by bringing 
each Vo m to agreement with V$ from the CF set. 
Thus we offset the OHP measurements by —4.97 



km s 1 , measurements from Boyajian et al. ( 2005 1 



by +2.7, and measurements from McSwain et al. 
p007a| by +2.4. Note that the unaltered OHP 



http: / /physics. nist.gov/PhysRefData/ASD/index. html 



IDL is a registered trademark of Research Systems, Inc. 



measurements are presented in Table [T] We then 
repeated the period search and orbital fit with the 
corrected V r , which improved our orbital solution 
significantly. Figure [2] shows the results from the 
final period searches. Note that for the CLEAN 
method, peaks in the signal indicate likely frequen- 
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cies, while for the PDM method the minima indi- 
cate more probable frequencies. Our systematic 
V r correction removed f a in the CLEAN search, 
but it remains present in the PDM periodogram. 

We present the final orbital solution in Table 
[2] and the corresponding radial velocity curve in 
Figure [3j The mean spectrum of HD 15137 from 
our CF and OHP runs, shifted to its rest wave- 
length according to our orbital solution, is shown 
in Figure [4] The final orbital solution of HD 
15137 indicates a highly eccentric (e = 0.62) bi- 
nary with a low velocity semiamplitude (K = 13.6 
km s _1 ). These orbital elements are very similar 



to our preliminary orbit (McSwain et al. 2007a I, 



although the period is quite different. While the 
09.5 III primary is likely a massive star, the very 
low mass function of this binary suggests a low 
mass companion. Based on the runaway nature of 
the binary, the eccentric orbit, and the low mass 



function, Boyajian et al. (2005) proposed a super 



nova ejection scenario and a NS companion in HD 
15137, even though the system was not a known 
X-ray source. They proposed that HD 15137 may 
be a "quiet" HMXB, too widely separated for the 
NS to accrete a significant mass of stellar winds 
to produce the bright X-ray flux commonly as- 
sociated with X-ray binaries ( Liu et aT7||2006 1 . To 
investigate this quiet HMXB scenario, we describe 
XMM-Newton observations of the system below in 
§4- 

3. Physical Parameters from Spectral 
Models 



In McSwain et al. (2007a I, we measured the 



stellar parameters of HD 15137 using the Tlusty 
OSTAR 2002 grid of line-blanketed, non-LTE, 
plane-parallel, hydrostatic atmosphere model 



spectra for O-type stars (Lanz & Hubeny 20031 



The OSTAR2002 grid uses a microturbulcnt ve- 
locity Vf = 10 km s , and we assumed solar 
abundances for the star. In our earlier paper we 
found Vsini = 234 km s~\ T cff = 29700 K, and 
logg = 3.5 for HD 15137 using fits of the He II 
line profiles. Given that we have obtained higher 
resolution blue spectra for a new comparison, we 
repeated the spectral modeling using the Tlusty 
OSTAR2002 grid QLanz fc HubenypObT l. 

We measured V sin i by comparing the ob- 
served He I AA4143, 4388, 4471, 4713 line pro- 



files in the mean CF and OHP spectra to the 
OSTAR2002 model line profiles convolved with 
a limb-darkened, rotational broadening function 
and a Gaussian instrumental broadening func- 
tion. We determined the best fit over a grid of 
Usini values by minimizing the mean square of 
the deviations, rms 2 . The formal error, A Tasini, 
is the offset from the best-fit value that increases 
the rms 2 by 2.7rms 2 /iV, where N is the number of 
wavelength points within the fit region. We mea- 
sured a new value of Vs'mi = 258 ± 20 km s _1 , 
consistent with our earlier result. 

We measured T e g by comparing the observed 
He ILHe I equivalent width ratios to the equiva- 
lent width ratios in the broadened model spectra 
(iWalborn & Fitzpatrick||1990|). The OSTAR2002 



grid overlaps with the Tlusty BSTAR2006 grid 
over the range of T e g and log g appropriate for 
a late O-type star, although the BSTAR2006 grid 
uses a lower V t (2 km s _1 for our values of logg; 
Lanz fc Hubeny|2007 1 . To consider whether lower 



values of V t may be appropriate for HD 15137, we 
also compared our spectra to both grids over the 
appropriate range in T e fj and log g. We measured 
equivalent widths, W\, using numerical integra- 
tions over each line profile, and we estimate an 
error of 10% in each observed W\ due to errors 
in the continuum placement and intrinsic noise. 
By fixing logg at each of three possible values 
(logg = 3.25, 3.5, and 3.75), we determined a 
nominal value of T c ff from our observed W\ ratios 
by interpolation with the model W\ ratios. The 
resulting values of T c g and log g are listed in Ta- 
ble [3] Using the mean T e g for each value of log g, 
we then compared the broadened model spectrum 
to our mean observed spectra over the full wave- 
length range and determined a reduced % 2 from 
the difference. The resulting values of reduced 
X 2 for each fit are listed in Table [3] Based upon 
our renewed inspection of the He I and He II line 
strengths, our new measurements for T G g — 29700 
K and logg = 3.5 agree perfectly with our earlier 
study. We estimate the error in T e g to be 1700 
K based on the standard deviation of the three 
measured values, and we estimate an error in logg 
to be 0.25 dex due to the spacing of the grids of 
model spectra. 

The final broadened model spectrum from the 
OSTAR2002 grid is presented in Figure |4| Com- 
pared to the model spectrum, HD 15137 has a 
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slightly weak C III A4187 line and strong N III 
AA4379,4514 lines. This may be evidence that the 
atmosphere is enriched with CNO-processed gas. 
Finally, we note that there is no evidence for an op- 
tical companion present in our observations based 
upon the excellent agreement of our observations 
and this model of an isolated O-type star. 

4. XMM-Newton Observations 

We observed HD 15137 with the XMM-Newton 
observatory on 2008 August 3, observation ID 
0553810201, for approximately 20 ks. Based on 
our proposed orbital solution, these observations 
took place at orbital phase cf> = 0.69. The field 
of view was centered on HD 15137 for a direct on- 
axis view. The three EPIC cameras were operated 
in Full Frame mode with the medium optical fil- 
ter, and the Optical Monitor was turned off during 
the observation due to the brightness of the star. 
The RGS instruments were operated in Standard 
Spectroscopy mode, but due to the faintness of the 
source the RGS data were not used. 

The EPIC Observation Data Files (ODFs) and 
event lists were provided by the standard XMM 
Pipeline Processing System. Using the XMM- 
Science Analysis System (SAS) version 7.1.2, we 
filtered the event lists using the standard cutoff 
of 0.35 counts s _1 for the MOS cameras and 5 
counts s _1 for the pn camera using the evselect 
command to exclude times of high particle back- 
ground. The resulting effective exposure times 
(good time intervals) for each camera are listed 
in Table H 

We extracted the source spectra from the event 
lists using a circular region with radius of 25" for 
the MOS cameras and 32" for the pn camera. To 
extract the background spectrum for each camera, 
we used a partial annulus region with inner and 
outer radii of 80" and 110", respectively, centered 
on the source. Because another possible weak X- 
ray source was detected in this annular region, 
we excluded the portion of the annulus near that 
source. The resulting count rates for the source 
and background regions are listed in Table g] As 
a weak X-ray source, we found no indication of 
pile-up in the observation. 

The resulting EPIC spectra of HD 15137 have 
low signal-to-noise (S/N), but they are consistent 
with a soft thermal source typical of isolated O- 



type stars (Sana et al. 20061. We fit the two 



MOS spectra simultaneously, over the range 0.5- 
2.3 keV, using a variety of warm absorbed (wabs; 
Morrison & McCammon 1983 1, single temperature 



(1-T) and two temperature (2-T) thermal mod- 
els available with Xspec version 11.3.2ag, includ- 



ing Raymond-Smith (Raymond & Smith 19771, 



MEK AL ( |Mewe et al.||1985[ |1986| |Liedahl et al. 
1995[), APEC (|Smith et al.|2001|), bremsstra hlung 



(Karzas & Latter 1961 Kellogg et al. 1975), and 



plain blackbody models. We repeated our fits for 
the pn spectrum using the same group of models. 
A sample 1-T fit of the pn spectrum is shown in 
Figure [5] The resulting fits were equally good for 
the 1-T and 2-T models, but a statistical F-test re- 
veals that the second temperature component does 
not significantly improve the fits. We also cannot 
distinguish between the quality of the various 1-T 
thermal model fits due to the low S/N. We weakly 
constrain the temperature to 0.10 < kT < 0.25 
keV, and the neutral hydrogen column density to 
2.8 x 10 21 <nH < 8.6 x 10 21 atoms cm" 2 . There 
is no evidence of any hard X-ray photons detected 
in our data. 

Our measured value of nH for HD 15137 can 
be compared to its reddening, determined from 
the ultraviolet, optical, and near-infrared spectral 
energy distribution (McSwain et al. 2007a I. In 
that work, we measured reddening E(B — V) = 
0.43 and a ratio of total-to-selective extinction, 
R = 3.18. The optical extinction is thus Ay = 
R x E(B — V). Several authors have published 



nH - A 


v relations 1 


Bohlin et al. 


1978 Wolk et 


al. 


2006 


Giiver & Ozel 


2009 


1 that predict nH = 



2.5 — 3.2xlO zi atoms cm -2 based on the measured 
reddening of HD 15137. Our measured range in 
nH is consistent with Ay, although our warm ab- 
sorbed 1-T models also allow values of nH that 
exceed the expected value by as much as a factor 
of 3. 

Unfortunately, interpreting our X-ray spectra 
is made somewhat ambiguous since the emergent 
spectra of cooling NSs is soft and nearly ther- 
mal, very similar to a blackbody (Treves et al. 



2000). We performed fits of the EPIC spectra 



using a model of the hydrogen atmosphere of a 



NS (nsa; |Pavlov et al.||1991| |Zavlin et al.||1996 |. 

As recommended by those authors, we fixed the 
NS mass to Mjvs = 1-4 Mq and the radius to 
Rns — 10 km. We used a fixed magnetic field 
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strength fixed to B — 10 12 G and included warm 
absorption fixed to the predicted nH — 3 x 10 21 
atoms cm~ 2 based on the observed Ay. Fitting 
the two MOS spectra simultaneously results in an 
unredshifted effective temperature of the NS of 
T cffjNS = 437000±lf ° ° Q ° K and a normalization of 
1.72 x 1CP 5 . From the fit of the pn spectrum, we 
obtain T offiNS = 955000^^^™ K and a normal- 
ization of 7.08 x 10~ 8 . The reduced x 2 values are 
equally good as the 1-T models described above, 
although the NS temperatures from the MOS and 
pn spectral fits are very inconsistent with each 
other. 

In order to place an upper limit on any hard 
power law component that may originate from an 
accreting compact companion, we repeated the 
warm absorbed, 1-T models with an additional 
power law component with photon index T = 2 
(fixed). All of the best fit parameters from the 1- 
T fits, including nH, fcT, and their normalizations, 
were fixed in Xspec. We then refit each model, al- 
lowing only the normalization of the power law 
component to vary. In every case, the best fit nor- 
malization for the power law component was zero. 
We then used the steppar routine to investigate 
the 90% confidence limit for the power law nor- 
malization. Upon removing the thermal compo- 
nent from the models, we used the fixed T and the 
upper limit for its normalization in the remaining 
absorbed power law model to determine the up- 
per limit for the X-ray flux, Fx, of the putative 
compact object. Fits of the pn and MOS spectra 
indicated unabsorbed Fx < 10~ 14 erg cm~ 2 s _1 
(over the range 0.2-10 keV). At a distance d = 2.42 
kpc ( McSwain et al.|2007a ), this corresponds to an 
X-ray luminosity Lx < 10 31 erg s _1 for any hard 
power law component. 

5. Discussion 

We can loosely constrain the mass of the com- 
panion in HD 15137 using the mass function, 
/(m), from our new orbital fit. Assuming the 
09.5 III star has a mass of 20.6 M Q (Martins et al. 



2005|, and noting that there is an 87% probability 
that the binary inclination i > 30°, we constrain 
the companion mass to 1.4 < M 2 < 3.0M Q . A 
plot of allowable values of Mi and M 2 is shown 
in Figure [6j The most probable range in M 2 cer- 
tainly suggests that a low mass NS or BH could 



present in the system. The kick velocity from a 
prior supernova event is a reasonable explanation 
for the current high peculiar space velocity (63±12 
km s _1 ; McSwain et al.|[2007a l for the system. A 
white dwarf companion is unlikely due to the ini- 
tial mass of the supernova progenitor. 

HD 15137 has a m easured mass loss rate , 
logM = -6.5 M yr~ x ( |Howarth fc Prinja||l989 1 
and a terminal wind velocity, Voo = 1690 km s _i 
( |McSwain et al.||2007b[ ). Using these wind param- 
eters and our new orbital solution, we can pre- 
dict the X-ray luminosity emitted from a possible 
wind- accreting NS in the system using a simple 



Bondi-Hoylc wind accretion model (Lamcrs et al 



1976). With the O star mass and a NS compan- 
ion with mass 1.4 M Q , we should expect to ob- 
serve an unabsorbed Fx ~ 2 x 10~ 12 erg cm 2 s _1 
(corresponding to an unabsorbed X-ray luminos- 
ity logL x = 33.2 erg s _1 ) at (f> = 0.7 if wind 
accretion was indeed present. A more massive, 
accreting BH should result in a proportionately 
higher Lx- However, the XMM-Newton obser- 
vations rule out such a bright X-ray source by 
at least two orders of magnitude, clearly ruling 
out an accreting BH companion. O-type stars 
typically exhibit order-of-magnitudc variations in 



their mass-loss rates (McSwain et al. 2004), so if 



the mass loss rate has decreased with time since 
measured by iHowarth & Prinjal ((1989), the ex- 



pected wind accretion luminosity should decrease 
proportionately. However, variability in M is not 
sufficient to explain the discrepancy between the 
predicted Fx and our observed upper limit. 

The dynamical age of the putative NS in HD 
15137 also allows us to discriminate against the 
cooling NS model. A young NS cools to a tem- 
perature ~ 10 6 K after about 10 3 -10 4 yr after its 
birth, and it is expected to remain at that roughly 
constant temperature for ~ 10 5 yr (Treves et al. 



2000). Since the binary was ejected 10 yr ago 



from the open cluster of its birth (Boyajian et al. 
|2005[ ), then it should have cooled to temperatures 
far below our measured T e g NS. 

If there is a NS or BH companion present in 
HD 15137, it is more likely in a quiescent state. 
Quiescent NSs in the propellor regime have been 
observed with hard X-ray spectra with T ~ 2 and 
unabsor bed luminosities of 33.5 < log Lx < 37.5 



erg s (Campana et al. 20021. Our upper limits 



on the hard X-ray flux of HD 15137 clearly rule 
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out such a NS in the propellor regime. NSs in the 
ejector regime are potentially observed as radio 
pulsars, but our pulsar search did not detect any 



such candidate (McSwain et al. 2007b I. We can- 



not firmly rule out a NS companion in HD 15137, 
but our observations suggest that a companion is 
unlikely to be present. However, a quiescent BH 
is quite possible. Black hole X-ray binaries spend 
most of their lives in quiescence with an X-ray 
luminosity as little as 10 -6 of an active accretor 
( |Pszota et al.||2008||Cor"iat et al.|2009| , consistent 
with our upper limits on Lx- 

If HD 15137 has a "twin" companion with mass 
ratio q = M 2 /Mi ~ 0.95 ( |Pinsonneault fc Stanek 
2006; Lucy] 2006 1, the system must have an im 
probably low i ~ 6°. The late O-type or early B- 
type companion would have spectral lines heavily 
blended with the primary star, making it difficult 
to detect from V r variations but possibly identi- 
fiable using doppler tomography. Boyaj ian et al.| 
( 2005 ) failed to detect an optical companion in HD 
15137 using doppler tomography but we repeated 
such an effort here with the advantage of a new or- 
bital solution for the binary. We assumed a wide 
range of 0.05 < q < 0.95 in an effort to detect a 
massive or low mass secondary. The reconstructed 
spectra from doppler tomography do a poor job 
of matching realistic H B aimer line profiles, and 
all possible reconstructed spectra maintain con- 
stant He I:He II line strength ratios. These line 
ratios would be expected to differ from the pri- 
mary star in any cooler optical companion, so we 
cannot claim any detection of a possible optical 
secondary star with the tomography. 



Boyajian et al. (2005) noted that the travel time 



since the ejection of HD 15137 presents an un- 
usual paradox since it is longer than the expected 
lifetime of an O-type star. Perets (20091 offers 



a possible explanation: the close binary may have 
been ejected from its parent cluster by a dynamical 
ejection, while later mass transfer rejuvenated the 
massive star and extended its lifetime. The orbital 
period and probable high mass ratio of HD 15137 
are consistent with the strong mass transfer sce- 
nario ( |Perets| [2009 ) . However, with a total binary 
mass almost certainly > 20M Q , the binary might 
have originally been composed of two ~ IOM0 
stars. Perets predicts a very small parameter 
space for such massive binaries that will experi- 
ence strong mass transfer; most of the massive bi- 



naries in their simulations experience weak mass 
transfer or mergers. While rejuvenation through 
mass transfer is possible, it seems an unlikely ex- 
planation for the current properties of HD 15137. 
On the other hand, it is also possible that HD 
15137 experienced a supernova after the binary's 
ejection via dynamical processes. In such case, 
the age of the compact companion would not cor- 
respond to the dynamical age found by |Boyajian| 
et al. (20051, and a young, cooling NS compan- 
ion might be present. However, such a two-step 
evolutionary scenario would obscure the true dy- 
namical history of HD 15137; any supernova kick 
velocity would alter the system's trajectory and 
render any measurement of its travel time unre- 
liable. It is impossible to draw any conclusions 
regarding the dynamical ejection and subsequent 
evolutionary history of HD 15137. 

6. Summary and Conclusions 

We have searched for a compact companion 
in the massive runaway binary HD 15137, to no 
avail. XMM-Newton EPIC spectra of the system 
are consistent with a soft thermal source typical of 



isolated O-type stars (Sana et al. 2006). The dy- 
namical age of HD 15137 implies the system is too 
old to contain a cooling NS with T c ff,NS ~ 10 5 — 10 6 
K, although our X-ray spectra are consistent with 
such a NS. Our lack of detection of hard X-ray 
photons is inconsistent with the expected emis- 
sion from a quiescent NS in the propellor regime 
or a stellar wind-accreting NS, but a quiescent BH 
remains a possibility. 

Similarly, we have been unable to identify an 
optical secondary in the system. Doppler tomog- 
raphy based on our updated orbital solution was 
unsuccessful, and spectral modeling of the primary 
is consistent with a single 09.5 III star. Weak N 
enrichment and C depletion are observed in the op- 
tical spectra, possible evidence of CNO processed 
gas at the stellar surface. The companion proba- 
bly has a mass in the range 1-3 Mq, implying a 
very low mass ratio of q < 0.2. For a MS compan- 
ion, this would correspond to a late B star which 
would contribute a negligible flux compared to the 
bright O-type star. Such an optical companion 
would be difficult if not impossible to detect in 
the close, eccentric orbit of HD 15137. 

Although it might be tempting to rule out a 
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compact companion based on the lack of its de- 
tection in our XMM-Newton observation, it is dif- 
ficult to imagine how the binary might have been 
ejected as a runaway with a low mass stellar com- 
panion in a close gravitational encounter. HD 
15137 bears many similarities with the class of 
Be/X-ray biniaries (BeXRB): a moderately long 
orbital period, high eccentricity, low mass func- 
tion, and a rapidly rotating primary star. Al- 



though Boyajian et al. (20051 reported no Ha 



emission from HD 15137, the rapid rotation and 
possible nonradial pulsations make the star a good 
candidate for a transient Be star that experi- 



ences sporadic disk disappearances (McSwain et 
al.|2009[ ). Many BeXRB are X-ray transients that 



arc invisible in X-rays except at times when the 
NS encounters higher density disk gas (there are 
no known BH companions in BeXRB systems; |Bel-| 
czynski & Ziolkowski 2009 ) . X-ray outbursts may 
happen either when a Be star disk outburst oc- 
curs and the disk density increases close to the 
NS, or when the NS crosses the plane of the disk 
(in cases where the orbital and disk planes are dif- 
ferent). The XMM-Newton observations may have 
occurred when conditions were unfavorable for ac- 
cretion. 
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Table 1 

Radial Velocity Measurements of HD 15137 



HJD V r <j 

(-2,400,000) Telescope (km s" 1 ) (km s" 1 ) 



53652.642 ... 


OHP 


-39.4 


3.8 


53654.538 ... 


OHP 


-37.5 


2.4 


53980.585 ... 


OHP 


-55.7 


2.6 


53982.552 ... 


OHP 


-42.7 


2.2 


53984.545 ... 


OHP 


-43.5 


1.8 


54033.406 ... 


OHP 


-53.6 


9.3 


54033.495 ... 


OHP 


-56.9 


10.3 


54033.602 ... 


OHP 


-54.4 


8.8 


54034.670 ... 


OHP 


-56.5 


9.2 


54035.413 ... 


OHP 


-58.0 


1.3 


54035.437 . . . 


OHP 


-59.8 


4.9 


54035.458 ... 


OHP 


-56.9 


8.8 


54035.480 ... 


OHP 


-59.8 


6.1 


54035.502 ... 


OHP 


-54.4 


10.0 


54035.524 . . . 


OHP 


-57.2 


6.7 


54035.545 ... 


OHP 


-59.6 


2.8 


54035.568 ... 


OHP 


-59.6 


3.2 


54035.588 ... 


OHP 


-56.2 


5.2 


54035.611 ... 


OHP 


-53.5 


4.4 


54035.632 ... 


OHP 


-56.2 


5.2 


54035.655 ... 


OHP 


-52.4 


6.3 


54396.439 ... 


OHP 


-29.1 


5.9 


54397.516 ... 


OHP 


-31.5 


3.5 


54400.646 ... 


OHP 


-31.7 


4.6 


54401.598 ... 


OHP 


-31.6 


5.9 


54402.588 ... 


OHP 


-28.4 


5.2 


54405.599 ... 


OHP 


-25.4 


5.2 


54406.590 ... 


OHP 


-22.8 


2.6 


54407.532 ... 


OHP 


-33.6 


1.7 


54408.591 ... 


OHP 


-29.1 


4.0 


54409.527 . . . 


OHP 


-29.7 


4.2 


54410.475 ... 


OHP 


-23.5 


2.7 


54411.551 ... 


OHP 


-27.2 


8.0 


54412.504 . . . 


OHP 


-28.6 


7.9 


54413.584 . . . 


OHP 


-35.9 


6.0 


54414.530 ... 


OHP 


-31.1 


7.0 


54415.545 ... 


OHP 


-34.8 


7.5 


54416.453 ... 


OHP 


-30.3 


4.6 


54417.467 . . . 


OHP 


-32.6 


6.1 


54418.383 ... 


OHP 


-33.3 


6.8 
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Table 1 — Continued 



HJD V r a 

(-2,400,000) Telescope (km s" 1 ) (km s" 1 ) 



54419.412 . . . 


OHP 


-45.3 


5.4 


54421.525 . . . 


OHP 


-54.5 


1.1 


54422.420 . . . 


OHP 


-56.4 


1.8 


54423.322 . . . 


OHP 


-51.0 


3.6 


54756.820 . . . 


CF 


-73.8 


17.9 


54757.753 . . . 


CF 


-57.3 


7.7 


54758.831 . . . 


CF 


-52.0 


4.1 


54759.771 . . . 


CF 


-51.1 


4.2 


54760.781 . . . 


CF 


-50.5 


2.1 


54761.815 . . . 


CF 


-50.6 


3.2 


54762.773 . . . 


CF 


-39.8 


4.7 


54763.779 . . . 


CF 


-48.0 


5.8 


54764.761 . . . 


CF 


-48.0 


2.1 


54765.741 . . . 


CF 


-45.2 


1.1 


54766.730 . . . 


CF 


-46.0 


4.5 


54767.955 . . . 


CF 


-47.3 


5.3 


54768.697 . . . 


CF 


-49.0 


11.2 


54769.721 . . . 


CF 


-48.4 


4.0 


54771.663 . . . 


CF 


-39.4 


1.4 


54771.674 . . . 


CF 


-41.1 


6.5 


54771.685 . . . 


CF 


-36.7 


1.4 


54771.696 . . . 


CF 


-37.4 


3.5 


54771.712 . . . 


CF 


-40.3 


5.2 


54771.722 . . . 


CF 


-36.4 


4.8 


54771.733 ... 


CF 


-39.5 


5.3 


54771.744 . . . 


CF 


-39.6 


2.7 


54771.754 . . . 


CF 


-38.1 


2.9 


54771.765 ... 


CF 


-34.4 


5.3 


54771.776 ... 


CF 


-38.9 


5.0 


54771.786 ... 


CF 


-40.8 


2.7 


54771.797 . . . 


CF 


-43.3 


2.9 


54771.808 ... 


CF 


-37.3 


6.7 


54771.818 ... 


CF 


-43.3 


4.9 


54771.829 ... 


CF 


-45.8 


4.2 


54772.641 ... 


CF 


-40.7 


7.5 


54774.751 ... 


CF 


-35.9 


4.8 


54775.702 ... 


CF 


-45.9 


3.8 


54777.696 ... 


CF 


-38.4 


11.4 


54779.754 . . . 


CF 


-43.9 


2.4 


54780.886 ... 


CF 


-32.8 


10.5 
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Table 1 — Continued 



HJD 






a 


(-2,400,000) 


Telescope 


(km s" 1 ) 


(km s- 1 ) 


54781.654 . . . 


CF 


-36.5 


11.6 


54782.629 ... 


CF 


-35.4 


4.2 


54783.659 ... 


CF 


-43.1 


5.0 


54784.654 . . . 


CF 


-39.3 


2.9 


54785.630 ... 


CF 


-43.0 


5.3 


54786.600 ... 


CF 


-47.8 


4.5 


54787.660 ... 


CF 


-44.3 


6.2 


54788.695 ... 


CF 


-42.0 


6.1 


54789.907 . . . 


CF 


-35.8 


2.2 


54790.909 ... 


CF 


-34.0 


10.8 


54791.943 ... 


CF 


-32.0 


10.3 



Table 2 

Orbital elements of HD 15137 



Element Value 



P (d) 55.3957 ± 0.0038 

T (HJD - 2,400,000) 54421.991 ± 0.064 

if(kms^) 13.56 ±0.15 

lofkms- 1 ) -42.06 ±0.11 

e 0.6239 ± 0.0088 

lu (°) 152.17 ±0.86 

f(M 2 ) (M Q ) 0.00685 ± 0.00030 

aisini(i? Q ) 11.60 ±0.17 

r.m.s. (km s" 1 ) 4.28 
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Table 3 

He Line Ratios and Effective Temperature 





, W A (HeII) 
iU & W x (Uel) 


Toff (K) 

I lug (J — O.ZO 1 


Tcff (K) 
i lug (j — o.ou ) 


Tcff (K) 

l lug (J — o. to) 


Using OSTAR2002 grid: 










He II A4542:He I A4471 


-0.656 


27041 


27814 


28719 


He II A4542:He I A4388 


-0.295 


28695 


30153 


30971 


He II A4200:He I A4144 


-0.122 


30099 


31202 


32390 


Mean T cS 




28600 


29700 


30700 


Reduced x 2 from CF spectrum . 




11.0 


2.1 


9.6 


Reduced x 2 from OHP spectrum 




4.3 


1.4 


7.9 


Using BSTAR2006 grid: 










He II A4542:He I A4471 


-0.656 


27764 


28889 


29954 


He II A4542:He I A4388 


-0.295 


29983 


31527 


33532 


He II A4200:He I A4144 


-0.122 


31640 


33798 


36645 


Mean T cff 




29800 


31400 


33400 


Reduced x 2 from CF spectrum . 




2.8 






Reduced x 2 from OHP spectrum 




1.7 







Table 4 

Journal of XMM-Newton observations of HD 15137 





Time of Mid Exposure 


Performed 


Effective 


Source Count 


Background Count 


Camera 


(,ID-2,400,000) 


Duration (ks) 


Duration (ks) 


Rate (cnt s _1 ) 


Rate (cnt s _1 ) 


MOS1 


54681.56704 


19.667 


15.90 


0.0100 ± 0.0008 


0.0036 ± 0.0008 


MOS2 


54681.56741 


19.672 


17.50 


0.0069 ± 0.0006 


0.0014 ± 0.0007 


pn 


54681.57333 


18.032 


18.10 


0.134 ± 0.003 


0.018 ± 0.003 
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0.8 



4465 4470 4475 
WAVELENGTH (A) 



Fig. 1. — Variations in the He I A4471 line pro- 
file of HD 15137 during the CF long stare. The 
spectra have been smoothed with a Gaussian of 
FWHM — 5 pixels and are offset vertically for clar- 
ity. 
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Fig. 3. — Radial velocity curve of HD 15137. 



The single point from Conti et al. (19771 is plot- 
ted as an open square, points from Boyajian et 



(2005) as open triangles, points from McSwain 



et al. (2007a I as open diamonds, points from the 



OHP in this work as filled circles, and points from 
the CF in this work as filled squares. Our long 
stares were performed at <j> w 0.0 and <fi w 0.3. A 
typical error bar, assuming a — 7 km s _1 in addi- 
tion to the intrinsic V r error of ±5 km s _1 due to 
rapid line profile variations, is also shown. 



0.04 



Fig. 2. — Frequency search results using PDM 
(top) and CLEAN (bottom) methods. The sam- 
pling frequency / s , alias frequency f a , and true sig- 
nal frequency / are also marked in the top panel. 
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Fig. 5. — XMM-Newton pn spectrum of HD 
15137. The spectrum has been fit in XSPEC us- 
ing a warmly absorbed thermal model, wabs(apec ), 
with nH = 6.3 x 10 21 atoms cm~ 2 , kT = 0.235 
keV, and a normalization of 2.083 x 10~ 4 . The 
reduced x 2 = 1-33 for this fit. 
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Fig. 4. — Mean spectrum of HD 15137 from our 
CF observations (A < 4450 A) and from OHP (A > 
4450 A). The CF spectrum has been smoothed 
for presentation. A model spectrum interpolated 



from the Tlusty OSTAR2002 grid ( |Lanz fc Hubeny 
2003), with Vsini = 258 km s"\ T eff = 29700 K, 
and log g = 3.5, is shown for comparison. 
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Fig. 6. — Mass diagram of HD 15137 is plotted 
for a range of inclination angles (solid lines) . The 
expected Mi from Martins et al. ( 2005| ) is also 



shown (vertical dashed line) . 
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